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ABSTRACT 



Swift J164449+573451 is a peculiar outburst which is most likely powered by 
the tidal disruption of a star by a massive black hole. Within the tidal disruption 
scenario, we show that the periastron distance is considerably smaller than the 
disruption radius and the outflow should be launched mainly via magnetic activ- 
ities (e.g., Blandford-Znajek process) otherwise the observed long-lasting X-ray 
afterglow emission satisfying the relation L x oc M can not be reproduced, where 
Lx is the X-ray luminosity and M is the accretion rate. We also suggest that 
Lx oc M may hold in the quick decline phase of Gamma-ray Bursts. 

Subject headings: accretion, accretion disks-black hole physics-Gamma Rays: general- 
radiation mechanism: non-thermal 



1. INTRODUCTION 



Swift J164449.3 +573451 (Sw J1644+57) triggered the Swift B urst Alert Telescope (BAT) 
on 2011 March 28 (jCummings et al. 1 1201 ll : iBurrows et al.l 120111 ). As revealed by the late 
optical observations, this transient lay at the center of a galaxy at the redshift z = 0.3534 
( jLevan et al. ll2011al Jb[h In the first few weeks, the average isotropic luminosity in the 0.3- 
10 keV band was about 10 47 — 10 48 ergs _1 . Several months later, it was still about a few 
xl0 45 ergs _1 , well above the Eddington limit. The X-ray emission declined as Lx oc t~ 5 ^ 3 
during the time interval from 10 5 s to 10 6 s after the tr ig ger, while the subsequent decline can 



be a pproximated as L x oc t" 4/3 flLevan et al. IboilaU bl: 
201 it ). 



Bloom et al. 2011 



Cannizzo et al. 
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The super-long duratio n (> 8 months) o f the X-ray activities essentially rules out the 
Gamma-ray Burst models (jShao et al. Il201ll ). Inst ead it strongly favors the model of tidal 



disruption of a (gian t ) star by a massive black hole (IBloqm et a 



Levan et al. 2011b; Cannizzo et al. 2011; Shao et al 



2011 



2011 



Burrows et al. 2011 



Wang fc Chend 12012 ). As 



for the central black hole (BH), it is impossible to measure the mass directly. Indirect 
constraints indicate that the mass should be in the ranee 1O 6 M - 10 7 M© (detailed analysis 
can be seen in Cannizzo et al. 2011). In the tidal disruption scenario, if a star passes within 
the disruption radius Rrp ~ R*(M BH j M*) 1 / 3 , the BH's tidal gravity exceeds the star's self- 
gravity and consequently the star is disrupted, where i?* (M*) is the radius (mass) of the 
disrupted star and Mbh is the mass of the BH. After about hours to weeks, part of the 
remnants remains on bound and will return to the perice nter of the orbit where the material 



start to be accreted inwar d, releasing a flare of energy (IRees I Il988l . Il990t iPhinney 
Strubbe fc Quataert II2OO9L 



1989 



Theoretical calculations of the tidal disruption events suggest that the immediately 



accreted unbin ding gas falls back to the pericenter at a rate Mfb oc t 5 / 3 (IRees I Il988 



Phinnev 1989), and the subsequ ent disk accretion follows a rate M^oct 4 / 3 if the disk is 



thick (jCannizzo fc Gehrels II2009I ). So for t > 10 5 s, the X-ray emission luminosity of Sw 



J1644+57 is proportional to the accretion rate, i.e., Lx oc Mfb- In this work we pay special 
attention on the physical implication of such a relation. The physical parameters of the 
disrupted star are also investigated. 



2. The physical parameters of the disrupted star 



There is an express, derived in the Newto nian limit, f or the timescal e of return of the 
most bound stellar material to the pericenter (IRees Ill988l ; IPhinney Ill989l ) 



t 



fb 



2nGM BH (2AE) 



-3/2 



0.048 yr( 



M 



BH 



5 x 10 6 M 



lV2, 







■ M* 



R Q 



3/2^3 



(1) 



where AE = k^ JL {^ff L ) 1 / z , the dimensionless coefficient k depends on the spin-up state 
of the star. If the star is spun up to the break-up spin angular vel o city, we have k 
3. If the spin- up effect is negligible then we have k ~ 1 (IRees 1 1 19881 ; I Aval et al. 1 1200 
The dimensionless coefficient /i = R p /Rt is taken to be a free parameter in the following 
discussion, where R p is the periastron distance of the star. 

The star could be spun-up via tidal interaction. In linear perturbation theory, the 
spin- up angular velocity is given by (|Press fc Teukolsky 1 119771 ; lAlexander fc Kumar 1 12001 
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Alexander & Livio Il200ll : iLi et al. 1 120021 ) 



u a ^ T 2 (/x 3 / 2 ) 
u p 2Ifi 3 



(2) 



where u p = v p /r p is the orbit angular velocity of the star at the pericenter, / is the stellar 
momentum of inertia in units of M*Rl, T 2 is the second tidal coupling coefficient and depends 
on the structure of the star and the eccentricity of the orbit. For an n = 1.5 polytrope 
star of mass 0.76 M & and radius 0.75 R Q , Alexander & Kumar(2001) found / ~ 0.21 and 
T 2 (l) ~ 0.36, corresponding to u) s /u p ~ 0.86 for /x = 1. Furthermore, they showed that 
the numerical simulations including nonlinear effects led to a larger energy transfer from the 
orbit to the star and a larger spin-up than that predicted by linear theory. Therefore we 
take k = 3 as the fiducial value in our analysis. For completeness we also present the results 
in the non-spinning case (i.e., k = 1). 



The bound material returns to pericenter at the rate (jPhinney II1989I ) 
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2AM ,t - 1. 



fb 



-5/3 
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where t s is the time of initial tidal disruption, AM is the mass that falls back to pericenter 
and the dimensionless factor / is defined as / = AM/M*. 

When the accretion rate of the fall-back material is highly super-Eddington, only a 
fraction (1 — f out ) of such material forms a disc and can be accreted all the way down 
to the central BH, i.e., M — (1 — jout)M^. The remaining part will i nstead leave the 
system undergoing a strong radiation pressure. IStrubbe fc Quataert I (120091 ) took a constant 
fout— 0-1- However, numerical simulation indicates that the para meter f mit is a growin g 
function of M fh /M Edd} reaching f mt ^0.7 for M fb /M Edd = 20 JDotan fc Shavivl boioh . 
Considering the observed X-ray luminosity Lx ~ 10 47 — 4 x 10 48 erg s _1 during the first 10 6 
s after the trigger (Burrow et al. 2011), even for a radiation efficiency as high as 0.1, we 
need an accretion rate M = lOLx/c 2 ~ 5 x 10~ 7 — 2 x 10 _5 M Q . The Eddington luminosity 
can be scaled as L Edd ~ 6.25 x 10 44 M 6 7 erg s _1 , so the Eddington rate M Edd = 10L Edd /c 2 = 
3.5 x 1O _7 M s _1 . We then have M/M Edd ~ 1 — 60. We define a free dimensionless parameter 
£ = f(l — fout), the fraction of the material that is actually accreted onto the central BH. 

Assuming that the jet radiation efficiency is e during the stellar debris fallback accretion, 
the intrinsic jet luminosity then can be given by 



eMc 2 w 2.5 x 10 42 erg s - 1 (^)(^)( 



M 



BH 



1 1/3, 



0.0r v 0.05 yv 5 x 10 6 Mt 
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where the efficiency is normalized to e ~ 0.01 (In sec 3.3 we will show the jet should be 
launched mainly via magnetic activities (e.g. B-Z effect) and the efficiency is about 0.01). 
As shown below, the conclusion drawn in this section is independent of the value of e. 

When most of bound debris falls back to the pericenter, the jet luminosity peaks at 
t — t s = tfb and can be estimated as 

L jpeak « 3.9 x 10 44 erg s'^— )(JL)( Mb » )-i/2 A2 A-3/2 - 3 (5) 
J,p & V 0.01 A 0.05 A 5 x 1O 6 M Mq Rq P v; 

The observed maximal X-ray luminosity is ~ 4 x 10 48 erg s _1 . For a collimated emitting 
region with a half-opening angle #j, we have the constraint Lx, P cak > 2 x 10 46 (^-) 2 erg s -1 . 
Assuming that most of the radiated energy is in the X-ray band during the time interval, we 
have 

W j W l 0.05 K 5x 106M Q j WW ^ ~ W 

The observed X-ray fluence Sx (0 s < t < 10 7 s) suggests a total energy Ai?x = 
9? j Lx{t)dt/2 ~ 1 x 10 51 (p J j-) 2 erg. The total mass of the accreted material is thus 

The approximated mass-radius relationship can be scaled as (7^) = (^)' ? - F° r the main 
seque nce stars, we have n 0.8 for O.1M < M# < 1M and 77 « 0.6 for 1M < M* < 
lOM (jKippenhahn fc Weigert Ill994l ). With equations (6) and (7) we obtain 

V ') 1 - 3 " / V 3 ( - M **- )' 1/2 >46.4, 



'M ' ^ v 5xlO 6 M 



which then yields 



^ < 0,36( Mb} ! )-i/«(j^L)V3-W2. (9) 
P ~ v 5 x 1O 6 M V M v ; 

Interestingly, the parameter /1 is independent of e, £ and Its dependence on both the stellar 
mass and the black hole mass is also very weak. The above analysis is under the condition 
that the star is spun-up to the break-up spin angular velocity, i.e. k = 3. For the non- 
spinning case k = 1, we can obtain a more stringent result \i < 0.16( 5x A ] / ^6f j - Q )^ 1 ^ 6 (^) 1 ^ 3- ^ 2 - 
Therefore, we conclude that the periastron distance is likely well within the t idal disruption 
radius (i.e., it is a plunging event), in agreement with ICannizzo et al. I ( l201ll ) . 

Based on the observation data, the peak accretion rate can be derived with Eq.5 and 
Eq.7, 

Mpeak = = eM : Lx ' peak = 4.4 x 1O" 5 M g-V J_)(_^) (10) 

peak ec 2 AE X V O.O5 A 1.1M V ; 
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The accretion rate M can be scaled as M 
process. 



peak \ 



■ t-u 

lyr 



5//3 during the fall-back accretion 



In this plunging event, the disrupted star's orbit is likely to mis-align with the equatorial 
plane of the spinning central BH. A tilted accretion d isk should be formed a nd the jet aligne d 



with the disk normal vecter is expected to precess (IStone fc Loeb 1 120121 ; iLei et al. II2012I ). 



Saxton et al.(2012) have analyzed the X-ray timing and spectral variability of Sw J1644+57 
and found the periodic modulation, possibly due to the jet precession. 



3. The radiation mechanism in the fallback phase 

After the time t ~ 10 5 s, the observed X-ray luminosity followed the fall-back accretion 
rate (Levan 2011a,b; Bloom et al. 2011), i.e. L oc M. Such a relationship have shed some 
light on the underlying physics. 



3.1. Thermal X-ray radiation from the disk? 



While the fallback accretion rate is super-Eddington, the stellar material returning to 
pericenter is so dense that it can not radiate and cool. The gas is most likely to form 
an advective dominated accretion flow (ADAF) accompanied with powerful outflow, which 
dominates the emission. Most of the radiation will be emitted from the outflow's photosphere. 
When the pho tosphere lies inside the outflow, the photosphere's rad ius and temperature can 
be written as ( jStrubbe fe Quataert 112009c iRossi fc Begelman 1120091 ) 



R 



pit 



^foutfv 



M 



'M. 



fb NR l/2 d 
- rt p,3H s - fl s 

Edd 



and 



T ph ~ i x io° k f :ru 



1/3 f 1/3, 



M 



fb 



Edd 



-5/12 M -l/4 D-7/24 
M 6.7 K p,3R s 



(12) 



where R s = 2 ™ Bg is the Schwarzschild radius, and f v is the ratio of terminal velocity of 
gas with the escape velocity at the radius ~ 2R P . However, the photons escape from the 
photosphere are mainly in the UV optical band and have a blackbody spectrum. The UV 
optical emission luminosity can be scaled as vL v ~ 4tt R 2 vh vB v {T vh ) oc R^ h T ph oc Mf b 19 ^ 2 . 

These optical photons could be Compton-scattered by the relativistic electrons in the 
outflow. The energy of the photons getting scattered is hvic ~ D 2 ^j 2 hv, where D = l/[r(l — 
f3cos8)} is the Doppler factor, T is the Lorentz factor of the outflow, 9 is the angle between 
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the outflow axis and the observer's line of sight, 7 is the Lorentz factor of the relativistic 
electrons. The energy of inverse Compton scattered photons can peak in the X-ray band if 
the parameter D ~ 1 and 7 ~ 10. However even in this case, the X-ray luminosity does not 
satisfy the relation L x oc M/&, inconsistent with the observational data. 



3.2. Neutrinos annihilation launched outflow? 



When the mass accretion rate is high enough, the accretion proceeds via neutrino cooling 
and neutrinos can carry away a significant amount of energy from the inner regions of the 
disk. The mechanism i s used to explain th e launch of at least some gamma-ray burst outflows 



(IPopham et al. I Il999l : Fan & Wei 



approximated by a simple formula ( jZalamea fc Beloborodov 1 1201 lh 



2011 



Fan et al. 



2012). The luminosity may be well 



K- v « Llxl0 52 X -f (f|£f 3/2 (13) 

L ign 

m g / 4 M ign < M < M trap } erg s"\ (14) 

where m = M/M Q s _1 , Xms = Rms{p) I Rs, Rms is the radius of the marginally stable orbit, 
M ign is the mass accretion rate to ignite the neutrino emitting, M trap is the mass accretion 
rate when the emitted neutrino becomes trapped in the disk and advected into the black 
hole. The characteristic accretion rates M ign and M trap depend on the viscosity parameter 



a and the mass of central BH. Based on the work of iBeloborodov I (120031 ). for //-transparent, 
the accretion rate should be as large as M ign > 7.6 x lO 30 ^) 1/2 (^)(^) 2 g s -1 . For this 
event Sw J1644+57, a plunging one, the peak accretion rate is M peak w 4.4 x 1O _5 M s _1 , 
which is far less than M ign . Hence we conclude that the observed X-ray emission could 
not be produced by neutrino annihilation (please see Shao et al. (2011) for an alternative 
argument disfavoring the neutrino mechanism). 



3.3. Poynting-flux dominated outflow? 



Extracting ener gy from the rotating black hole may be possible through the Blandford- 
Znajek mechanism (IBlandford fc Znajek 1119771 ). Such a process is based on the expectation 
that the differential rotation of the disk will amplify pre-existing magnetic fields until they 
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approach equipartition with the gas kinetic energy. For a black hole of mass Mbh and 
angular momentum J, with a magnetic field B± normal to the horizon at Rh, the power 
arising from BZ mechanism is given by (e.g. Thorne et al. 1986) 



8 An J ma . 



(15) 



where J„ 



= GM 2 /c is the maximal angular momentum of the black hole. The fac- 
tor Up = Qp(Qh — ^f)/^1 depends on the angular velocity of field lines Qp relative to 
that of the black hole, fit- Usua l ly we adopt ojf = 1/2, which maximizes the power 
output (IMacdonald fc Thorne I Il982l ; iThorne et al. I Il986l ) . We follow the common assump- 
tion that the magnetic field in the disk will rise to some fraction of its equipartition value 
P mn .g = |r~ ~ otP in the inner disk. The pressure P = pc 2 is given by ( Armitage fc Natarajan 



19991 ) 



V2M 



(5 + 2 £ ) 1 / 2 (GM) 1 / 2 r 5 / 2 , 



(16) 



Una 

where e is the parameter governing the property of the disk. For a thick disk we have e < 1 
otherwise e > 1. In the inner region of the disk, we assume B± B, Rk, R h = GM/c 2 . The 
BZ power for the case of a maximally rotating black hole (i.e., J = J max ) can be estimated 
as 

L BZ « 7 x 1(T 3 (5 + 2e) 1/2 Mc 2 1 (17) 

corresponding to an efficiency e^z — Lbz/Mc 2 ~ 10 -2 for the thick disk model. In the 
thin-disk scenario the radiative efficiency can be as high as ~ 0.1. For Sw J1644+57, at the 
time t — t s = tfb + 10 6 s, the mass accretion rate is about 3.3 x 1O~ 8 M s -1 , the observed 
luminosity is 2Lbz@J 2 ~ 5 x 10 47 erg s _1 (^j/0.1) -2 , consistent with the observation. 

Our conclusion that the outflow powering the super-long X-ray emission should be 
launched via magnetic activities (e.g., B-Z mechanism) is consistent with that of Shao et 
al. (2011). Lei & Zhang(2011) have also analyzed the jet launched by B-Z mechanism and 
then constrained the physical parameter of the central BH. One interesting finding is that 
the central BH should have a moderate to high spin. 

In the Poynting-flux domina ted outflow, the X-ray emissio n could be due to the dis- 
sipation of the magnetic field ( jUsov I Il994t iThompson I Il994l ). There are several mag- 
netic field dissipation models that could produc e the observed emission, such as the global 
MHD condition breakdown model ( jUsov Ill994j ). the gradual magne tic reconnection model 
( iDrenkhahn fc Spruit 1120021 ). the magnetiz ed internal shock mo del (IFan et al. 1 120041 ) . and 
the collision induced reconnection model ( IZhang &: Yan 1 1201 ll ). For illustration, here we 
take the global MHD condition breakdown model to calculate the emission. By comparing 
with the pair density (oc r -2 , r is the radial distance from the central source) and the density 
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required for co-rotation (oc r 1 beyond the light cylinder of the compact obje ct), one can 



estimate the radius at w hich the MHD condition breaks down, which reads (jUsov 1 11994 
Zhang fc Meszdros lliooi l 



Tmhd ~ 5 x 10 20 cm( 



L 

W 1 ' 



vl/2 



■W [ W 



)( 



10' 



(if 



where a is the ratio of the magnetic energy flux to the particle energy flux, T is the bulk 
Lorentz factor of the outflow, t v>m is the minimum variability timescale of the central engine. 
Beyond this radius, intense electromagnetic waves are generated and outflowing particles are 
accelerated (e.g. Usov 1994). Such a significant magnetic dissipation process converts the 
electromagnetic energy into radiation. 



At tmhd, the corresp onding synchrotron radiation frequency can be estimated as (jFan et al. 



2005; Gao k Fan 200 



V m ,MHD 
13(2-2)1 



1.5 x 10 18 Hz 



1 + f 'plf £ )f 2 ( ° ) 3 ( ^ \f tv,rn - 

i.35 ' l o.r pl io J 4o n io 2 ' 



(19) 



where C p = (^)[ 3(p _ 1} 



e e is the fraction of the dissipated comoving magnetic field energy 
converted to to the comoving kinetic energy of the electrons, and the accelerated electrons 
distribute as a single power-law dri/d^e cx 7~ p , £ < 1 reflects the efficiency of magnetic 
energy dissipation. So most energy is radiated in the X-ray band. 



4. Clue to the X-ray steep decline following the prompt emission in 

gamma-ray bursts 

The observations of Sw 1644+57 suggest that the long-lasting fall-back accretion onto 
a black hole can produce energetic X-ray emission and the radiation luminosity traces the 
accretion rate (i.e., Lx cx M). One interesting question is whether similar process takes place 
in Gamma-ray Bursts (GRBs) or not. The answer may be positive. Here we just discuss 
the X-ray steep decline (quicker than t~ 3 , see Fig.l of Zhang et al. (2006) for illustration) 
following the prompt emission in GRBs. In the collapsar model, a fraction of the gas in the 
core of the collapsing star has not sufficient centrifugal support and directly forms a central 
black hole. The rest of the material will have sufficient angular momentum to go into orbit 
around the black hole. The fall-back accretion rate is tightly related with the pre-collapse 
stellar density profile, which is of the form p cx r~ T . Numerical simulation has found when 
the outermost 0.5 M & layer of the star (wher e r > 5) is accreted , the fall-back accretion rat e 



can get to the value M/f, cx t 3 or steeper (IMacFadyen et al. I l200ll ; iKumar et al. 1120081 ). 
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If the relation Lx oc M still holds, one has an X-ray emission decline steeper than t 3 , in 
agreement with the observational data0. 

5. Conclusion and Discussion 

Swift J164449+573451 is a peculiar outburst which is most-likely powered by the tidal 
disruption of a star by a massive black hole. In this work we find out that the ratio of the 
periastron distance to the disruption radius \i < 1, implying that SW J1644+77 is a plunging 
event (see section 2). The mass of the plunging star however can not be tightly constrained 
due to its strong dependence on the poorly understood parameters e (the radiation efficiency), 
£ (the fraction mass of the star that is eventually accreted into the central) and 9j (the half 
opening angle of the collimated outflow). 

As a tidal disruption event, the accretion rate M at late times (say, t > 10 day) is 
relatively well understood and is widely believed to be oc t~ 4//3 . The detected X-ray emission 
Lx shows a rather similar decline behavior. Since the forward shock origin of the long- 
lasting and highly variable X-ray emission has already been convincingly ruled out (Shao 
et al. 2011), the X-ray emission has to be from an outflow launched by the accreting black 
hole. These two facts strongly suggest that Lx oc M, which can shed valuable light on 
the underlying physics, in particular the energy extraction process. Three kinds of possible 
mechanisms have been examined and only the Poynting-flux dominated outflow model is 
found to be able to account for the data (see section 3 for details). Therefore the magnetic 
activity at the central engine (e.g., Blandford-Znajek process or B landf or d- Payne process) 
plays the main role in extracting the rotation energy of the black hole and then launching the 
outflow. We suggest that L x oc M may also hold in the quick decline phase of Gamma-ray 
Bursts. 
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